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Abstract — Hydrodynamic journal bearings are critical 
power transmission components that are carrying increasingly 
high loads because of the increasing power density in various 
machines. Therefore, knowing the true operating conditions of 
hydrodynamic journal bearings is essential for machine design. 
Oil film thickness is one of the key operating parameters 
describing the operating conditions in hydrodynamic journal 
bearings. Measuring the oil film thickness in bearings has been 
a demanding task and therefore the subject has been studied 
mainly by mathematical means. 

The aim of this study was to determine the oil film thickness 
in real hydrodynamic journal bearings under realistic operating 
conditions. The study focused on engine bearings. Calculations 
were carried out to determine the oil film thickness and to 
understand its relationship with other operating parameters. 
Analytical calculations are highly complex and very difficult to 
perform using MS Excel. So Mathcad is a good solution for 
performing this type of complex calculations and formulations. 

Once the template is ready, it can be used for other engine 
bearings to determine the minimum oil film thickness by 
analytical means. After performing whole calculations, the 
template is used to integrate with actual CAD model of the 
components. The results can be used in the development and 
validation of mathematical methods for research into 
hydrodynamic journal bearings. Integration process for ball 
bearing is shown in detail using Mathcad and Creo. Then actual 
CAD model will be automatically modified as dimension in 
Mathcad worksheet modifies after successful regeneration in 
Creo parametric. Design template calculations results are in 
good agreement with results given by “KISSsoft” software. 

Index Terms — Hydrodynamic lubrication, Integration 
procedure between Mathcad and Creo, Journal bearing, 
Mathcad formulation, Oil film pressure, Oil film thickness. 


I. INTRODUCTION 

The power density in various machines, for example in 
internal combustion engines, is increasing year by year due to 
growing demands for mechanical and economic efficiency. In 
machine design, one of the consequences of an increase in 
power density is that critical power transmission components 
have to carry increasingly high loads. 
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Hydrodynamic journal bearings are typical critical power 
transmission components that carry high loads in different 
machines. In machine design, therefore, it is essential to know 
the true or expected operating conditions of the bearings. 
These operating conditions can be studied by mathematical 
means, for example in field or laboratory tests with engines 
and by calculation. 

Numerous studies of the operating conditions of 
hydrodynamic journal bearings have been made during the 
last decades. Still, the case is far from closed. For example, 
there are a limited number of studies that carry out an 
in-depth examination of the true operating conditions of 
bearings in true-scale experiments. There is also a need for 
experimental studies to verify the theoretical ones. 

The operating conditions of hydrodynamic journal 
bearings can be described by a set of tribological variables 
called key operating parameters. The key operating 
parameters most directly related to the bearing lubricant-shaft 
contact are the oil film temperature, oil film thickness and oil 
film pressure. These three key parameters can be determined 
by mathematical means with varying levels of complexity. 
Until now, oil film pressure in hydrodynamic journal bearings 
has been studied mainly by mathematical means, because the 
experimental determination of oil film pressure has been a 
demanding or even an unfeasible task. Under real operating 
conditions, there are typically many practicalities that 
complicate the experimental determination of true oil film 
pressure in a certain point or at a certain moment. The oil film 
may be extremely thin and therefore sensitive to different 
disturbing factors, for example defects in geometry. In 
addition, the level of the oil film pressure may be extremely 
high or have a high level of dynamic variability. 

II. BASICS OF THE OPERATION OF 
HYDRODYNAMIC JOURNAL BEARINGS 

Lubrication reduces friction between two surfaces (such as 
sliding surfaces of a bearing and a shaft) in relative motion. It 
is typically categorised as boundary, mixed and 
hydrodynamic lubrication, as shown in example by Heywood 
(1988), Becker (2004) [12] and Gleghorn and Bonassar 
(2008) [11]. When a journal bearing operates under 
boundary lubrication, the sliding surfaces of the bearing and 
shaft are practically in direct contact and friction is at its 
highest level. Lower friction levels are achieved through the 
use of mixed lubrication, where the sliding surfaces are 
partially separated by the lubricant, and of hydrodynamic 
lubrication, where the sliding surfaces are completely 
separated by the lubricant. 
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Duty parameter = Dynamic viscosity x angular speed / pressure 

Figure 1. Stribeck curve showing the friction coefficient as a 
function of the duty parameter under boundary, mixed and 
hydrodynamic lubrication. 

To illustrate how friction varies under different 
lubrication conditions, Stribeck curves have been used 
widely in different engineering sciences. In Stribeck curves, 
friction coefficient is presented as a function of a 
dimensionless parameter calculated from dynamic viscosity, 
angular speed and pressure (Figure 1). The above-mentioned 
parameter is typically called the duty parameter or Hersey 
number. The minimum of the friction coefficient is reached at 
the critical value of the duty parameter, at the dividing line 
between the mixed and hydrodynamic lubrication zones. 
Hey wood (1988) [12] presented a Stribeck curve for a 
journal bearing. Methods for the calculation of Stribeck 
curves were studied by de Kraker et al. (2007) [16]. They 
calculated the friction coefficient as a function of the journal 
frequency at different values of the projected bearing 
pressure. 


Housing 



Figure 2. Hydrodynamic journal bearing with the diameter D and 
the width B carrying the bearing load L. The shaft with the diameter d 
runs at the angular velocity oo s . An oil film in the lubrication gap 
separates the sliding surfaces of the bearing and shaft. The eccentricity 
e and the angle p determine the position of the point with the minimum 
oil film thickness ho. 

The bearing unit consists of the bearing, housing, shaft 
and supporting bearings. The housing in which the bearing is 
placed for testing typically has a simplified cubic or 
cylindrical design, but housings of real types have also been 
used. A large high precision roller bearing with both radial 
and axial load carrying capacity is a common supporting 
bearing type used in low and medium speed applications. 
These have been used, for example, by Savagkan et al. (2002) 
[3]. In high speed or extreme load applications, sliding 
bearings have been used as an alternative, for example by 
Okamoto et al. (2000) [3] and Zhou et al. (2004) [10]. 

A hydrodynamic journal bearing (Figure 2) is designed to 
operate normally under hydrodynamic lubrication, where 
hydrodynamic pressure (Figure 3) in the lubricant keeps the 


sliding surfaces of the bearing and shaft separated from each 
other. The hydrodynamic pressure is caused by the sliding 
motion. 


P P 



Figure 3. Distribution of the hydrodynamic pressure p in the oil 
film on the sliding surface of a hydrodynamic journal bearing with 
the diameter D and the width B carrying the bearing load L. The 
shaft with the diameter d runs at the angular velocity oo s . 

Hydrodynamic journal bearings are simple but critical 
components and numerous parameters influence their 
operation. Therefore, research into the operation of 
hydrodynamic journal bearings is typically a demanding task 
requiring extensive knowledge of machine design. 

III. DETERMINATION OF MINIMUM OIL FILM 
THICKNESS FOR CRANKSHAFT BEARING 

Knowledge for bearings can be obtained by performing 
field or laboratory tests on engines and by mathematical 
means. There are advantages and disadvantages to each of the 
above mentioned ways to collect information about bearings. 

The density p T either be measured or calculated by 
different standard methods. For example, Dzida and 
Prusakiewicz (2008) [17] measured densities and specific 
heat capacities of different oils. Effenberger (2000) [3] 
presented an equation for estimating specific heat capacity as 
a function of the density and temperature. By using the 
approximate values and equations presented by Affenzeller 
and Glaser (1996) [1], for typical engines oils with a density 
of about 875 kg/m3 and a specific heat capacity of about 1.8 
kJ/kgK, it can be calculated that the relative change in density 
is typically about -0.07 %/K and the relative change in 
specific heat capacity is typically about 0.2 %/K. In this study, 
it was assumed that the density varies linearly as a function of 
the temperature at a constant pressure, it was estimated that 
the relative change in density Ap/AT was -0.07 %/K, and the 
density p T was calculated by the following equation: 

Pr= Pref[ 1-0.0007 (r Tref )/K] . (1) 

Where, p REF = Reference Density in Kg/mm 3 . 

T = Service Temperature in K. 

Tref = Reference Temperature in K. 

K = Kelvin. 

The minimum oil film thickness was either simulated by 
the simulation software or calculated. The calculations are 
presented below. 

The calculation was made in three phases. In the first 
phase, the Sommerfeld number S 0 (a dimensionless 
parameter used in bearing performance calculations) was 
determined approximately by the following equation, based 
on the measurement data: 
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S 0 = FC 2 / DWt]co .(2) 

Where, F = Radial Bearing Load. 

C = Relative Bearing Clearance. 

D = Diameter of the Bearing. 

W = Width of the Bearing. 
r| = ji p = Dynamic Viscosity, 
p = Kinematic Viscosity, 
p = Density of Oil. 

co = Hydrodynamic Angular Velocity. 

The density p T in equation (1) was calculated as a function 
of the temperature. In the second phase, the relative 
eccentricity 8 was determined approximately as a function of 
the Sommerfeld number S 0 and the width-to-diameter ratio 
B/D of the bearing (see Appendix A). The approximation was 
made for a plain bearing with the width-to-diameter ratio B/D 
= 32 mm / 85 mm =0.376, and the relative eccentricity s was 
calculated by the following approximate equation for 1 < S 0 < 
200: 

£ = kX 2 .(3) 

Where, ki = Coefficient. 

S 0 = Sommerfeld Number. 
k 2 = Coefficient. 

The values of the coefficients k x and k 2 in the equation (3) 
are presented in Table 1. 

Table 1. Values of the coefficients ki and k 2 with different values of the 
Sommerfeld number S 0 . 



K 

k 2 

1 < So < 10 

0.798 

0.073 

10 < So <100 

0.897 

0.022 

100 < So <200 

0.980 

0.0028 


In the third phase, the following equation was used to 
calculate the minimum oil film thickness h 0 as a function of 
the bearing diameter, relative bearing clearance and relative 
eccentricity: 

K - [l/2]DC(l —f).(4) 

Where, D = Diameter of the Bearing. 

C = Relative Bearing Clearance, 
s = Relative Eccentricity. 

Due to the use of a simple calculation method, it can be 
estimated that the relative error in the calculated minimum oil 
thickness was high, about +10%, when the results are 
compared to values determined by detailed calculation 
methods. 

The oil used for measurement of minimum oil film 
thickness is SAE 15W-40 with following lubricant 
properties: 


Table 2. Properties of SAE 15W-40 oil. 


Properties 

Reference 

Specification 

Viscosity grade 

SAE J 300 

15W-40 

Density at 20 °C 

ASTM D 1298 

875 kg/m 3 

Viscosity at 100 °C 

ASTM D 445 

14.5 mm 2 /s 

Viscosity at 40 °C 

ASTM D 445 

110 mm 2 / s 

Pour point 

ASTM D 97 

-27 °C / -17 °F 

Flash point 

ASTM D 92 

224 °C / 435 °F 


The procedure for determination of minimum oil film 
thickness by using PTC Mathcad Prime software is as 
follows: 
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THEREFORE. PERCENTAGE ERROR IN BOTH RESULT IS 3.2328% 


Figure 4. Determination of minimum oil film thickness using 
PTC Mathcad Prime software. 
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IV. DESIGN METHODOLOGY FOR MATHCAD 

In this section we will discuss the design methodology 
used to create the procedure for the diesel engine components 
using Mathcad software. The design process for 
determination of minimum oil film thickness for plain journal 
bearing consists of following steps: 

Consider, 

F = Radial load on bearing in N. 
c = Clearance in bearing in mm. 

D = Diameter of bearing in mm. 
d = Diameter of journal in mm. 

C = Relative clearance in bearing. 

W = Width of bearing in mm. 

Tref= Reference temperature in K. 

Pref = Density of oil at reference temperature in 

kg/mm 3 . 

T m = Oil inlet temperature in K. 

Tout = Oil outlet temperature in K. 

T = Service temperature in K. 

p T = Density of oil at service temperature in 

kg/mm 3 . 

p = Kinematic viscosity in mm 2 /s. 
r| = Dynamic viscosity in Pa-s. 
co = Speed of revolution in rpm. 

S 0 = Sommerfeld number. 

s = Relative eccentricity. 

h 0 = Minimum oil film thickness in pm. 

Ki, K 2 = Constants as per the range of 
Sommerfeld number. 

Diameter of journal is given by - 
d = (D -c) 


Due to the use of a simple calculation method, it can be 
estimated that the relative error in the calculated minimum oil 
thickness was less, about +5%, when the results are compared 
to values determined by detail calculation methods. 

V. MATHCAD TO CREO INTEGRATION 

Mathcad can be integrated with Creo Parametric. There is 
provision for integration in Creo Parametric in Analysis tab 
located in command bar of software. This section explains 
how to integrate Mathcad with Creo Parametric. Following 
process steps gives brief idea for integration of Mathcad to 
Creo Parametric by using crankshaft ball bearing, because it 
is difficult to simulate oil film thickness in journal bearing 
using Creo. 

Create a ball bearing defined by outer and inner diameter 
with specified width in Creo Parametric. 


-POO 



TP• H a 


Figure 5. Ball bearing with 50 mm inner diameter. 

Go to tools tab and click parameters tab. Add parameter 
as scale to 1 and click ok. 


Relative clearance in bearing is given by - 
C — d d 

Service temperature is given by - 

T ~ {Tin + T out ) / 2 

Density of oil at service temperature is given by - 

Pt = PrefC 1 - 0 - 0007 ! 7 ’-^)^] 

Constants for 15W-40 oil at standard working condition 
[17] are given by - 

B = jlog ^log ^ / [mm / 5 )+0.7 j j - log j^log {// 2 / [mm 2 1 s j+0.7 j j J / {log [T 2 I K} - log (7| 

A = {log [log / (mm 2 / s) + 10.7)] j + B {k>g (T 2 / A")} 
Kinematic viscosity of oil is given by - 

r r 10 (A- B [io g (r/^)])i 1 

jli = \ 10 -0.7 \mm Is 

Dynamic viscosity of oil is given by - 
r| = pp T 

Sommerfeld number is given by - 

S 0 = FC 2 / DWr/co 

Relative eccentricity is determined by - 

s = ^Sq 2 

Minimum oil film thickness is given by - 
\ =[l/2]DC(l-f) 

The minimum oil film thickness obtained by Mathcad 
software is 3.2449 pm and the KISSsoft result is 3.14 pm. 
The percentage error between these two results is about 



Figure 6. Specifying scale in Creo. 


In analysis tab, go to prime analysis and then load 
Mathcad file for integration. 



3.2328 %. 
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After loading Mathcad file, assign scale to 1 and the 
dimension which you want to change by integration. Here, we 
use bearing diameter to 48mm. Assign these parameters as 
input and output parameters. Then save the Mathcad 
worksheet. 
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Now prime analysis tab shows mapped value and then 
click compute. Add feature name for analysis and close it. 


DETERMINATION OF MINIMUM OIL FILM THICKNESS FOR 
PLAIN JOURNAL BEARING. 


SCALE -1 vwuvowwta 

RADIAL LOAD ON BEARING ^ •- 

CLEARANCE IN BEARING _ 

DIAMETER OF BEARING 
DIAMETER OF JOURNAL 
RELATIVE CLEARANCE IN BEARING 


SCALE= 1 
Ft- 25000 N 
ct= 0.055 mm 
D-SCALE -48 mm 
d-.=(D-c) =0.0479 n 
Cs=-=0.00115 


Figure 8. Mathcad worksheet with 48 mm inner diameter. 

Now go to Creo prime analysis interface. Tab shows the 
name of Mathcad file loaded in it. Click Auto-map and create 
Creo parameter to prime mapping. Right click and Select 
Creo parametric parameter. Parameter selection pop up and 
select scale parameter. 



Figure 9. Parameter selection and mapping. 
After this, right click and select prime input scale. 



Figure 10. Mapping prime input. 

Create prime to Creo parameter mapping. Right click and 
Select prime output parameter. Parameter selection pop up 
and select scale parameter as D. 



... & -“I s. 



Figure 12. Mapping results and computation. 

Go to tools tab and add local parameters in relations and 
assign relation to parameter for verifying relations. Specify 
relations in feature parameter. Click on local parameter. Add 
parameter to relations and verify relations. Click ok. 



Figure 13. Specifying relations and its verification. 

The software gives yellow indication to complete 
regeneration step. Then right click on extrude and sketch 
command applied and edit definition for checking 
dimensions. 


/ ca n. a e> x a & 



Figure 14. Regeneration of inner part to 48 mm diameter. 

Then complete regeneration and check dimension of inner 
part of bearing. It will shows new dimension which applied in 
Mathcad worksheet. Here, diameter of bearing is modified 
using integration process from 50mm to 48mm. 

Right click on ANALYSIS 1 command applied and edit 
definition with clicking next for changing dimensions of 
Mathcad worksheet. 

Click compute tab. The dimensions changes 
automatically. Regenerate the changing parameter step. Then 
check extrude and sketch dimension by editing definitions of 
both. 


■v an 


PO O 

. ^ 



Figure 11. Prime to Creo parameter mapping. 


Figure 15. Bearing with modified inner diameter to 48 mm. 
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VI. CONCLUSION 

The main result of this study was the mathematical 
formulation of realistic oil film thickness in real type 
hydrodynamic journal bearings at various operating points 
across the realistic operating range. This template can be used 
for other engine bearings, minimum oil film calculation. 

Analytical calculations are highly complex and very 
difficult to perform using MS Excel. So Mathcad is a good 
solution for performing this type of complex calculations and 
formulations. 

The minimum oil film thickness obtained by Mathcad 
template is 3.2449 pm and the KISSsoft result is 3.14 pm. 
The percentage error between these two results is about 
3.2328 %. 

After performing whole calculations, the template is 
used to integrate with actual CAD model of the components. 
In this paper, Integration process for ball bearing is shown in 
detail using Mathcad and Creo. Then actual CAD model will 
be automatically modified as dimension in Mathcad 
worksheet modifies after successful regeneration in Creo 
parametric. Mathcad template results are in good agreement 
with “KISSsoft” results. 

APPENDIX 

The relative eccentricity 8 was determined approximately 
as a function of the Sommerfeld number S G and the 
width-to-diameter ratio B / D. The approximation of relative 
eccentricity was based on a graph which is presented in a 
standard for hydrodynamic journal bearings with a static load 
(DIN 31652 Teil 2, 1983) and in which the Sommerfeld 
number S 0 is presented as a function of the width-to-diameter 
ratio B / D at different values of the relative eccentricity 8. 
The estimated relative eccentricity in a typical case (with the 
width B = 32 mm, the diameter D = 85 mm, and the 
width-to-diameter ratio B/D = 0.376)is presented in Figure 
16. 



Figure 16. Estimated relative eccentricity £ as a function of the 
Sommerfeld number S 0 . The width-to-diameter ratio B / D is 0.376, and 
1 < So < 200. 
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